Cooling of a small sample of Bose atoms with accidental degeneracy 
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A system of bosons in a harmonic trap is cooled via their interactions with a thermal reservoir. 
We derive the master equation that governs the evolution of the system and may describe diverse 
physical situations: laser cooling, symphatetic cooling, etc. We investigate in detail the dynamics 
of the gas in the Lamb-Dicke limit, whereby the size of the trap is small in comparison to the de 
Broglie wavelength of the reservoir quanta. In this case, the dynamics is characterized by two time 
scales. First, an equilibrium is reached on a fast time scale within the degenerated subspaces of the 
system. Then, an equilibrium between these subspaces is reached on a slow time scale. 
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I. INTRODUCTION 
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A. Bose— Einstein condensation 

The observation of effects related to the quantum sta- 
tistical properties of weakly interacting gases of atoms 
[jlj has become in the last decades a major challenge 
of atomic physics. Thus, a large part of the theoretical 
and experimental research has been focused during the 
90's on cooling atoms confined in traps at relatively high 
densities 0. These studies have led to the remarkable 
experimental realisation of a Bose-Einstein condensation 
(BEC) in Rubidium ||, and sodium vapors Evidence 
of BEC in a Lithium gas with attractive interactions 
has been also reported ||. These remarkable achieve- 
ments, have opened fascinating possibilities and applica- 
tions, such as the developement of a coherent source of 
atoms, or atom laser Q. 

The theoretical description of a system of ultracold 
bosonic or fermionic atoms is particularly convenient in 
the framework of second quantization. Quantum field 
theories of cold atoms [0, originally developed in a 
condensed-matter context, have been used in the diag- 
nostics of a Bose-Einstein condensate |§| , and in nonlin- 
ear atomic optics 

Several theoretical works in the recent years were more 
directly aimed to the dynamics of the cooling processes, 
dynamics of the possible phase transitions and the for- 
mation of quantum collective states. Those works con- 
cern both collisional cooling mechanisms, such as evap- 
oration Jl(],[ll| or symphatetic cooling Il2-M|, or laser 
cooling mechanisms, such as sideband cooling jl6), Ra- 
man cooling jL7|, and dark state cooling ]l8| . The latter 
processes allow to reach temperatures below the photon 
recoil energy E R [equal to (hk) 2 /(2M), where k is the 
laser wavevector and M the atomic mass]. One expects 
that the system might then enter a collective quantum 
state (such as Bose-Einstein condensate, or some ana- 
logue of it). In particular, under such conditions one 
hopes to realize also a coherent source of atoms. 



In general, the quantum dynamics of a system of cold 
atoms is a very complex many-body problem. Some of 
the above mentioned processes may be analyzed using 
quantum Boltzmann equations |]l9[ , po| . Starting from 
1994 a more general method based on a quantum mas- 
ter equation (QME) description J2l[] has been developed. 
In particular, a QME describing the dynamics of a small 
sample of laser cooled atoms in a harmonic microtrap has 
been proposed and analyzed |2^| . The quantum statis- 
tical nature of the atoms is reflected in the dynamics of 
the cooling process. In the case when the trapping poten- 
tials for the atoms in the ground and excited electronic 
states are different such QME might lead to multistabil- 
ity and generalized Bose-Einstein distributions p3[ | . The 
QME has the advantage that it permits to study atom 
number fluctuations in each of the trap levels, and thus 
provides a more complete description of the cooling pro- 
cess. In particular, it may be used, in principle, to de- 
scribe the dynamics of condensate formation. We have 
also derived a QME for symphatetic cooling |24|] and an- 
alyzed the possibility of achieving the condensation of a 
system of light particles which interact with a reservoir 
of heavier particles. In the context of nonlinear atomic 
optics a self-consistent Born-Markov-Hartree-Fock mas- 
ter equation has been derived and analyzed to study the 
spontaneous emission effects on atomic solitons pq| . 

Since that early works, the theory of QME for many- 
body systems has been very strongly developed. In par- 
ticular, C. Gardiner, P. Zoller and co-workers developed 
in a series of papers |^6| the QME and Quantum Ki- 
netic Boltzmann approach to describe the dynamics of 
the evaporative cooling. We have also extended the the- 
ory of the collective laser cooling to much more realistic 
situations, avoiding the microtrap assumption, and work- 
ing beyond the so-called Lamb-Dicke limit, in which the 
trap is of the size of the laser wavelength p7[ . In such 



1 



analysis, we avoided the reabsorption problem by work- 
ing in the so-called Festina Lente regime, in which the 
spontaneous emission rate is smaller than the trap fre- 
quency. 

One should stress that, apart from the area of quan- 
tum optics, master equations for quantum Bose or Fermi 
gases have been used in statistical physics p8[ . However, 
the master equation is usually postulated there starting 
from general statistical requirements. It describes the ap- 
proach towards the thermal equilibrium described by the 
Bose-Einstein or, correspondingly, Fermi-Dirac distribu- 
tions, it fullfills detailed balance conditions, and some- 
times it conserves some order parameters. The dynam- 
ics that it generates might have some universal proper- 
ties, but does not have a direct physical interpretation in 
terms of interaction with specific energy reservoirs. 

In most of the quantum optical examples, the mas- 
ter equation is derived via the elimination of the "bath" 
degrees of freedom starting from a more general theory 
that describes a very specyfic physical situation. The 
eliminated "bath" has a direct physical interpretation - 
it consists of photons, colliding atoms etc. Each of the 
jumps between the states of the system described by the 
QME usually corresponds to a very well defined physi- 
cal process of photon emission, absorption, atom-atom 
collision etc. p9[ . 

The QME approach, although valid in general, can 
only be used for practical calculations when the QME can 
be reduced to a set of kinetic equations, while the density 
matrix to a diagonal form. Such reduction is not always 
possible, and very often requires first the choice of an ap- 
propriate basis in the Hilbert space of the states of the 
system. In particular, the bare states of the system, that 
are the eigenstates of the system Hamiltonian in isolation 
from the "bath" , are not necessarily the right ones. For 
instance, in the process of cooling of an ideal bosonic gas 
confined in a microtrap, the reduction in the basis of bare 
states is possible only if additional assumptions are made. 
These statements seem a little surprising, since they hold 
even for arbitrary weak system- "bath" interactions, i.e. 
a situation in which the quantum Boltzmann equations 
should be valid ^0|. Note, however, that the validity of 
the quantum Boltzmann equations usually requires as- 
sumptions concerning quantum ergodicity, which simply 
do not hold in the above mentioned situations. The point 
here is that an ideal bosonic gas in a harmonic trap has 
plenty of degenerated bare energy levels. The quantum 
coherences between those levels (i.e. non-diagonal ele- 
ments of the density matrix) can survive for very long 
times, and contribute significantly to the dynamics. 



C. Degeneracy in many— body systems 

Let us enumerate by fh the eigenstates of a single atom 
Hamiltonian in the rotationally symmetric harmonic trap 
of frequency u>, where m is a natural number in one di- 



mension, a pair of natural numbers in 2D, a triple in 
3D etc. When we consider an ensemble of N atoms, the 
states of such an ideal gas can be written in the Fock 
representation as \n^,n^, . . .), where denote the oc- 
cupation numbers of the corresponding m-th eigenstatc. 
For noninteracting atoms there are two kinds of degen- 
eracies in such a system. First, a degeneracy of energy 
levels due to rotational invariance; that is for the states 
for which the sum of 's with a fixed sum of the compo- 
nents of m, is fixed itself. Obviously, such degeneracies 
are not present in ID. We shall not discuss them here. 
Second, there exists an accidental degeneracy, due to the 
particular symmetry of the harmonic potential p0[ | . This 
degeneracy occurs even in the case of ID: for instance for 
the states |0, 2, 0, . . .) and |1, 0, 1,0,...). Here, the state 
with two atoms in the first energy level has an energy 
2 x hio, which is equal to the energy of the state with one 
atom in the ground level and another atom in the second 
excited level (1 x 0hu> + 1 x 2Tiu>). Both kinds of degen- 
eracies are lifted up if one considers anisotropic trap with 
anharmonic energy levels. If one then assumes that the 
resulting energy level shifts are larger than cooling rates, 
one can evoke standard secular arguments to reduce the 
master equation to a diagonal form in the basis of the 
bare ideal gas states [|2|. It is precisely the subject of 
the present paper to study the situation in which such 
reduction is not possible, i.e. when the effects of acciden- 
tal degeneracy dominate the dynamics of the system. 

One could argue that such problem is purely academic 
since in real physical system atom-atom interaction will 
always lift the accidental degeneracy. The whole point 
is, however, that as long as the interactions arc not 
too strong, the system will still exhibit effects of quasi- 
degenracy. This will be the case when energy differ- 
ences between the quasi-degenerated levels will be small 
in comparison to Tiuj. Such condition is realised in not 
too dense systems, i.e. the system containing not too 
many atoms. The idealized theory presented in this pa- 
per is formulated for arbitrary number of atoms N, but 
in practice it applies to the situation when N is not too 
large (see Section VII). Nevertheless, in view of the com- 
plexity of the problem, it is in our opinion reasonable 
to treat the ideal case of non-interacting atoms in order 
to get insight into more realistic cases. That said, let 
us note at this point that in the last years the exter- 
nal modification of the s-wave scattering length (which 
dominates the atom-atom collisions at low energies) , has 
been theoretically investigated in different scenarios |u] , 
and also experimentaly demonstrated by employing the 
so-called Feshbach resonances [ f32"|]33"| . Remarkably, very 
recent experimental results show that a modification of 
the scattering length to very small values is experimen- 
tally feasible opening the fascinating possibility to 
acheive a quasi-ideal bosonic g that studied in the 
present paper. 
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D. Content of the paper 

The paper concerns thus the problem of cooling of an 
ideal Bose gas, or more precisely speaking a sample con- 
sisting of N atoms, in a perfectly harmonic microtrap. 
It should be noted that this problem is quite general. 
Atomic traps, although frequently anisotropic (see for 
example Ref. ||), can be designed to be, with a very 
good accuracy, harmonic. Moreover, even though in the 
small atomic samples atom-atom interactions will lift the 
exact degeneracy of energy levels, the system will remain 
quasi- degenerated. We expect that a cooled atomic sam- 
ple in such a harmonic microtrap will exhibit the effects 
of accidental quasi-degeneracy regardless of the method 
used for its cooling In order to stress this general as- 
pect of our study, we adopt here partially the statistical 
physics approach, and derive a master equation using a 
phenomenological model of the "bath" . In particular, 
our "bath" may represent one of the following two reser- 
voirs: Bl) a photon reservoir in the case of laser cooling 
pH; B2) an atomic reservoir consisting of other atoms 
in the case of symphatetic cooling p4|]. In the case B2) 
the bath atoms are trapped in a larger trap than the 
system atoms. Such situation (proposed and discussed 
in [24| ) could be realised if a small, say far-off-resonance 
dipole trap for system atoms was located inside a larger 
magnetic trap for the bath atoms Q]. 

Given that the resulting master equations for all these 
reservoirs have the same structure, the qualitative behav- 
ior for the cooling dynamics given by our specific model 
is quite general. The reason why we have chosen such 
a model for the bath is that it has the adventage that 
one can derive analytical formulas for the transition rates 
between different levels. We also stress that the mathe- 
matical treatment developed in this paper can be easily 
generalized to other physical situations, in particular for 
ultra-cold trapped polarised Fermi gases, which due to 
the suppression of the s-wave scattering induced by the 
Pauli principle, can be in an excelent approximation as 
ideal gases |g5| • 

This paper is organized as follows: In Section II we in- 
troduce the model, describing separately the system in a 
trap, the atomic bath, and the system-bath interactions. 
This Section is very much analogous to Section II of Ref. 
p4j , but is formulated for a different model of the bath. 
In Section III we derive the master equation governing 
the dynamics of the system, under Born and Markov ap- 
proximations. This equation is further analyzed in Sec- 
tion IV in the microtrap limit in terms of the Lamb-Dicke 
(LD) expansion, i.e. a systematic expansion in a small 
parameter n — ka, where a is the size of the trap ground 
state wavefunction, whereas A: is a typical wavevector of 
the bath quanta that is relevant for the cooling process. 
In Section V we discuss explicitely the breakdown of er- 
godicity due to the accidental degeneracy that occurs on 
a fast time scale in the lowest order of the LD expansion. 
Section VI is devoted to the discussion of the restoration 



of ergodicity on a slower time scale due to higher or- 
ders of LD expansion, whereas Section VII contains our 
conclusions. The paper also contains three appendices. 
Appendix A contains some useful formulas of the opera- 
tor algebra used in the paper, and describes the structure 
of the Hilbert space. Appendix B describes the details of 
the construction of some of the multiple vacua that ap- 
pear due to the accidental deneracy, whereas Appendix 
C presents matrix elements of the second quantized op- 
erators used in the paper. These elements are used for 
some calculations regarding the cooling rates. Finally, in 
Appendix D we present a proof of the decay of coherences 
on the slow time scale. 



II. DESCRIPTION OF THE MODEL 

We consider a system "A" of bosonic particles that are 
confined in a trap, and interact with a bath bosonic par- 
ticles "B". We assume that the particles "B" are prac- 
tically not affected by their interactions with the system 
"A" , so that "B" can be regarded as a phenomenological 
reservoir for "A" . The coupling to the bath represents 
the influence of some externally controlled cooling mech- 
anism (laser cooling, symphatetic cooling, etc.) on the 
system "A". The reservoir is assumed to be in thermal 
equilibrium at some given temperature. In this Section 
we introduce the Hamiltonian for the system "A" , the 
bath "B", and for their mutual interactions. The formal- 
ism is developed for the case of d spatial dimensions. 

A. Description of the system 

The system "A" is an ideal gas of N bosonic atoms 
of mass Ma confined by an isotropic harmonic poten- 
tial in d dimensions. In a second quantized form, and in 
the Fock representation the Hamiltonian describing the 
system can be written in the form 

Ha = J~] %uj{n x + n v + . . Qata^, (1) 
ft 

where u> is the trap frequency, n = (n x , n y , . . .) with 
n x , n y , . . . = 0,1,2,..., and at and a,i are creation and 
annihilation operators of particles in the n-th level of the 
harmonic potential, respectively. Since in the present pa- 
per we are interested in the effects of the accidental de- 
generacy, we neglect the contribution of the atom-atom 
interactions to the total Hamiltonian (for the discussion 
of its role see [p2|-p4[ , and Section VII) . 

B. Bath 

Similarly as in Ref. ^4j, we assume that the system 
"B" contains a practically infinite number of bosons of 
mass Mb embedded in a practically infinite volume, with 
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finite density ub- The free Hamiltonian for the bath "B" 
of particles in a second quantized form is 



1 



dke(k)b(k^b(k). 



(2) 



Here, A; is a wavevector in a d-dimensional space. The 
function e(k) represents the dispersion relation of the 
bath particles. For instance, for quasi-free particles it 
reads 



{hkf 

2Mb 



(3) 



where Mb is the effective mass of the bath "quanta" . In 
the following we shall use Eq. (||) but the theory is easily 
generalized to arbitrary shapes of the dispersion relation. 
In particular, we shall also consider massless bath quanta 
with a photonic-like linear dispersion relation. Note that 
such dispersion relation is in fact appropriate for both 
types of bath (Bl), and (B2) mentioned above. This is 
obvious in the case of laser cooling, but is also true for 
collisional cooling schemes, since no mass is produced or 
lost in the collision processes. The operators b{ky and 
b(k) are creation and annihilation operators of bath parti- 
cles corresponding to plane wave states with momentum 
k. In the case of laser cooling (Bl) k corresponds to a 
photon momentum, whereas in the case (B2) k is rather 
a momentum transfer associated with the system atom- 
bath atom collision. The operators b(ky and b(k) fulfill 
the usual commutation relations 



b(k),b(k') = b(ky,b(k 
6(fc),6(fc') t l =5 {d) {k-k')- 



0, 



In thermal equilibrium, the density operator describ- 
ing the state of the bath ps corresponds to the usual 
Bose-Einstein distribution (BED) |36|]. In this situation, 
we have 



(b(k)b(k')) = (6(fc)t&(fc')t) = o, 
{b{k)^b(k'))=n(k)5 {d) {k~k'), 



(4a) 
(4b) 



where n{k) is related to the number of particles with 
wavevector k, and is given by 



n{k) = 



-/3e(k) 



1 



-/3e(k) 



(5) 



In the above expression, [3 = 1/ (fcsT) is the inverse tem- 
perature, and z = exp(/3/u) is the fugacity, while /i de- 
notes the chemical potential. Note that /Lt = 0, z = 1 
for massless quanta, and for massive particles below the 
condensation point. Note also that both n(k) and e(k) 
only depend on k = \ k \ . Particle and energy densities are 
connected with these quantities by the relations 



n B = 



(27r) d 
1 

W) 



dkn(k), 
, i dkn(k)e(k), 



(6a) 
(6b) 



respectively. 



C. Interactions 



Within the present model the interactions between the 
particles and the bath describe the annihilation of an 
atom followed by its inmediate recreation, accompanied 
by absorption or emission of a single bath quantum. We 
have chosen such a model for the interactions since it 
is the simplest one that retains the effect produced by 
the accidental degeneracy in the cooling process. In any 
case, it may be regarded as a phenomenological inter- 
action that may be due to various physical mechanisms 
(system atom-bath atom collisions, laser-atom interac- 
tions, etc.). Similarly to the case of atom-atom collisions 
we employ here an analogue of the shape-independent 
approximation ]37| , |38| ] to write down the corresponding 
atom-bath quantum interactions. We assume that these 
interactions are local (i.e. have a zero range) in the spa- 
tial representation. Mathematically, this approximation 
means that the wave functions of both kinds of particles 
do not change significantly over distances characterizing 
interparticle potential in the relevant energy range. In 
the Fock representation the interaction Hamiltonian is 
given by 



Ha-b = ^2 J dk^nM'{k)a\a,n'b\k), 



where 



Ififi'ik) 



(2tt) 



dxipn(x)*'>Pn>(x) e ~ 



(7) 



(8) 



tpn(x) is the wavefunction corresponding to the n-th level 
of the harmonic oscillator, and k is a coupling constant. 
We have chosen the form of the interaction Hamiltonian 
([?]), because of its simplicity. We stress, however, that 
the qualitative (and to some extend quantitative) results 
of the paper do not depend on the particular choice of 
Ha-b- 

Without loss of generality, we can exclude from the 
integration over k in (Q) the value k = 0. This is clear 
since 



(o)4 a «' &t (°) 
(£4«*1 ( &t (°))' 



(9) 



is a constant shift operator of the zero momentum mode, 
proportional to the number of particles in the system 
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"A" . One can always perform a unitary shift transforma- 
tion of 6(0), and 6^(0) that cancels the term (0) and its 
hermitian conjugate in the Hamiltonian. Obviously, such 
transformation modifies the BED for the bath quanta 
with zero momentum, but the latter modification has no 
relevance for the transformed system-bath interactions 
in which the coupling to the bath zero mode is absent. 
Hence, from now on, in the integrals over k it will be im- 
plicitely assumed that k ^ 0. On the other hand, since in 
the next Section we are going to make a rotating wave ap- 
proximation (RWA) in the master equation derived from 
the Hamiltonian (0), we reduce Eq. (]?]) (as in Ref. 24 1) 
to the form 



Ha-b — H 



(H a +Hl). 



(10) 



Here, Ho contains the part of Ha-b given in (ffl) in which 
the sum is extended over values with J2s=x y T n s ~ n ' s ) = 
0. H a contains the part of Ha~b proportional to the 
bath annihilation operators b(k), with the sum over n, ft 1 
extended over the values for which ^2 S=X y (n s — n' s ) = 
a. 

The QME that we shall derive in the next Section 
will contain a Hamiltonian part (describing energy level 
shifts), and non-Hamiltonian part describing dissipative 
decay processes. It is worth mentioning that due to the 
RWA the Hamiltonian part of the master equation will 
not be generally correct |2l| , |39| . However, as in Ref. 
we will be only interested in the dissipative part of the 
master equation, which is correctly described under the 
mentioned RWA, provided the trap frequency lo is larger 
than the cooling rates. The latter assumption will be 
made all over the present paper. 



III. DERIVATION OF THE MASTER EQUATION 

The master equation for the above defined model can 
be derived following well-stablished procedures in the 
field of quantum optics pl] , p9[ , analogously to those dis- 
cussed in Ref. ]24jj . We first move to an interaction pic- 
ture defined by the unitary operator exp [— i{Ha + Hb)^\ ■ 
In this picture, the density operator p describing system- 
plus-bath degrees of freeedom fulfills the following equa- 
tion: 



dp(t) 
dt 



-[H A _ B (t),p(t)], 



(11) 



where the tilde indicates that the operators are ex- 
pressed in the interaction picture. We integrate for- 
mally this equation, and substitute back into (11). 
Next, we define the reduced density operator for sys- 
tem "A", pa — Ttb(p), where Ttb stands for the trace 
over the bath states, and make use of the fact that 
Ti b{[H a- b (0), p(0)]} = 0, since we assume that the 
density operator for the bath /Cb(0) is diagonal in the 



Fock basis (with respect to Hb), whereas Ha-b does not 
contain any diagonal matrix elements [the reader should 

recall that we have extracted the terms with k = in 
Ha-b, see Eq. (0)]. We obtain the following equation: 

dp~A{t) _ 
dt 

1 f 1 

/ dTTl B {[HA-B{t),[HA-B(t-T),p(t-T)}]}. (12) 

n Jo 

In the following step we perform Born-Markov approx- 
imation. For this approximation we have to assume that 
the correlation time r c of the reservoir is much shorter 
than the typical time over which p A (t) changes, i.e. the 
cooling time fj9f . The Born-Markov approximation is 
also related to the fact that the bath quanta are practi- 
cally not affected by their interactions with the system; 
this allows us to write p(t — r) = p A (t — r) <E> Ps(0). 
From the technical point of view, the correlation time t c 
can be defined as a time for which the integrand of ( p^ ) 
practically vanishes. For specific physical models of the 
bath it can be directly evaluated (see, for instance, Ref. 
|p4j|). The cooling time, on the other hand, depends on 
the physics of the interactions between atoms and bath 
quanta. It is controllable, and thus is assumed to be the 
longest time scale of the problem. In this case, we can 
safely substitute pA(t — r) by p~ A {t) m the integral Jl2]), 
and extend the upper limit of the integral to infinity. 

In the next steps we make use of Eq. ( |l0| ) and perform 
the RWA, i.e. neglect terms rotating at multiples of the 
trap frequency. Again, this approximation is based on 
the assumption that trap frequencies are large in com- 
parison to the cooling rates. 

Finally, taking into account the bath properties (Q), 
and coming back to the Schrodinger picture we obtain 
the following master equation: 



dpA 
dt 



■[H A + H' A _ A ,p A ] + Cp A , 



(13) 



where H' A _ A is the Hamiltonian term produced by the 
elimination of the bath in the master equation. Physi- 
cally, this term accounts for the energy level shifts due 
to the effective interaction between system particles via 
their interactions with the bath quanta. H' A _ A may 
also include the original atom-atom interactions (pro- 
vided they were present in the original Hamiltonian) in 
the spirit of independent rates approximation | |3S| ] - Simi- 
larly as in the case of the atom-atom collisions, all of the 
shifts caused by H' AA may be neglected in some situa- 
tions depending on the specyfic model of the bath, the 
size of the trap, and the number of atoms in the system 
(for details see [p2|-p4|). We shall omit them in the fol- 
lowing, and come back to the discussion of their role in 
Section VII. We shall therefore postulate a QME restrict- 
ing our attention to the Liouvillian C that describes the 
cooling process, 
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Q#0 



(14) 



where 



C-aRA — 



y . 

/ j n,n 
n.n'rh.rh 



- alaffata^pA - p A ata^/ata A /). (15) 

The sum in this expression is extended to n, n',m,m' 
fulfilling 



E 

s—x,y. 



(n s - n' s ) = a, 



E 

s—x : y.. 



(m s — m' s ) = —a. (16) 



Liouvillian ( |l5| ) acounts for transitions of particles from 
one level of the harmonic oscillator to another, experienc- 
ing a change in the energy of ahuj, and a corresponding 
absorption or emission of a bath quantum. Thus, the 
term with a = that conserves the energy does not en- 



ter ( 14) since we have excluded the bath quanta with zero 
energy from the interaction (see Sec. II). The terms with 
a > (a < 0), on the other hand, describe processes 
increasing (decreasing) the energy. These transitions are 
characterized by 



(17) 



x[n(k) + l]S[e(k) ~ huja], 
for a positive, and by 



xn(k)S[e(k) — fnj\a\ 



(18) 



for a negative. 

We stress here the fact that for any kind of interactions 
between the atoms and the bath, the corresponding Li- 
ouvillian has the same form as in (14) and (|l5|), provided 
in each interaction act one atom is annihilated and an- 
other created. This is the reason why our results of the 
following sections can be extended to other kind of inter- 
actions. Note, however, that in such cases the coefficients 
( |l7| ) and ( |l8| ) may have a much more complicated form 
(see, for example, Ref. 24 1). 



IV. LAMB-DICKE EXPANSION IN ID 

In this Section we perform the Lamb-Dickc (LD) ex- 
pansion of the master equation (pT[). From now on we 
shall concentrate on the one dimensional case and skip 
the vector notation. 

The LD expansion is valid in the situation when the 
bath quanta relevant for the cooling process have mo- 
menta k much smaller than the inverse of the size of 



the trap, a. Their corresponding wave functions [oc 
exp(±ikx)] vary slowly on a scale of a, and can be then 
expanded in Taylor series around x ~ 0. Since, accord- 
ing to Eqs. (|l7j ) and ([l8|) the relevant bath quanta have 
energies e(k) = huia, their corresponding momenta may 
be determined from the dispersion relation. For example, 
for the case of massive free particles characterized by the 
dispersion relation (||), the validity of LD at relatively 
low temperatures requires that 



2M B u\a\ 



1/2 



2M a uj 



1/2 



( M B \a\ 
\ M A 



1/2 



< 1, 

(19) 



i.e. the effective mass of the bath particles must be much 
smaller than that of the system atoms f4Q] . It is easy to 
find analogous conditions corresponding to other forms 
of the dispertion relation, including the case of massless 
bath quanta. In case of the laser cooling (Bl) at low tem- 
peratures the condition is simply that ka = 2-Ka/X < 1, 
where k (A) is the laser photon momentum (wavelength) . 

In the case of symphatetic cooling (B2) the condition 
is Aka < 1, where Ak is the typical momentum trans- 
fer in a collision act. In such collision the kinetic en- 
ergy of the incoming bath atom is of the order of 1//3, 
and its correponding momentum of the order of 2it/\b, 
where Xb denotes the thermal de Broglie wavelength. If 
the bath was cool, (3huj 3> 1, the momentum transfer 
would typically be equal to the final momentum of the 
bath atom. The validity of LD expansion would then 
require that the same condition as (1^) is fulfilled, with 
Mb denoting the real mass of the bath atoms, and that 
27ra/A B < (Mb/Ma) 1 / 2 < 1. 

In the LD regime we can expand the rates (|l7|), (118 
Denoting by 



k a — 



2M B uj\a\ 



1/2 



Va= k a a = y/ai], 
we obtain for n, n' , to, m' fulfilling Eq. (p 



(20) 
(21) 



2tt 

(k a ) 

x [n(fc Q ) + (l + sign(a))/2]-^ 

n Ken 



(22) 



and 

Tn,n' (k a ) 



2n 

,2 



Sn,n> -ifla (Vn+T^,n'-1 + V^n.n' + l) 

^ ( v / (« + 2)(n+l)J n ,„,_ 2 + (2n + l)5 n , n > 



+ y/n(n - l)8 n ,n'+2 ) + ■ 



(23) 



The master equation can be now rewritten using the 
above expressions in the interaction picture with respect 



G 



to Ha, and upon neglection of H' AA . It takes then the 
following form 



P = 



£ (0) + 4 X) + 4 X) + 0(oc T? 6 ) + ... p. (24) 



The "zeroth" order term is actually of the order of rj 2 
and has the form 



£ (0) p= T l [2A P A^ - A^Ap - pA^A] 
+r T [2A*pA - AA^p - pAA^] ; 



where 



r T = at 



r, 2 [n{ki) + 1] f 



(2tt) 2 ft 
(2^1 



v 2 KM] 



ft 2 fci 



whereas 



I 00 

^= -y= X! ^+l fl i fl n+l. 



1 V 

At= "y^ E + l4 +1 a„ 



(25) 

(26a) 
(26b) 

(27a) 
(27b) 



Note that quite generally: a) both cooling rates Tj, Tf 
are collective (i.e. proportional to TV); b) their ratio is 



n(k\ 



ze 



-phuj 



(28) 



Ti n(fci) + l 
Similarly, the higher order terms (of order oc 17 ) are 

£< X) p = -r x ry 2 [2 + CpA*) 

- (&Ap + A*Cp) - (p&A + pA^C)} 

- r T ?7 2 [2 (A*pC + C^pA) 

- (CA^p + AC^p) - ( P CA^ + P AC^)] , (29) 



with 



n=0 



G " = vm Eo+D 3/2 4 +1 « n , 

* Tl — 



(30a) 
(30b) 



and 



4 X) P= r[ 1} [2BpB t - B*Bp - pB^B] 



I 

r\ 1] [2B^pB - BB^p - pBB^} ; 



where 



,(i) 



N 



k 2 2ir f]\ 



vf = N 



(2tt) 2 ft 4 

k 2 2n r\\ 



n(fca) + 1] 



Mb 
ft 2 fc. 



(2tt) 2 ft 4 



(31) 

(32a) 
(32b) 



whereas 



^ oo 

B= 77^ E ^("+1)0 + 2)4^+2, 

^ Tl— 

1 X 

fit= -Tt? E VO + l)0 + 2)4 



+2 u n- 



(33a) 
(33b) 



n=0 



Note that similarly as in Eq. fl2q) the rates fulfill 



r T 1} _ n(k 2 ) = ze -20hu 



r W n(fc 2 ) + l 



(34) 



Finally, Eq. ( pif ) contains terms of higher orders oc rf ,rf , 
etc. 



V. ACCIDENTAL DEGENERACY AND THE 
BREAKDOWN OF ERGODICITY 

In this Section we discuss the dynamics governed by 
the lowest order term in the Lamb-Dicke expansion, i.e. 
by the equation 



P = 

First we observe that 



£ (o: 



p. 



[A,At] = l, 



(35) 



(36) 



i.e. the operators A, and A' represent an abstract har- 
monic oscillator Weyl-Heisenberg algebra, whereas the 
Liouville-von Neuman superoperator ( p5[ ) describes in- 
teraction of that harmonic oscillator with an effective 
thermal bath (see, for instance 21 1). The inverse tem- 
perature of the bath is given by 



hujf3 e 



log 



= Tilo/3 — (3p. 



(37) 



Note that the effective temperature is never greater than 
the temperature of the bath, (3 e > (3, since p < 0. Al- 
ternatively, one may view the above equation, as if the 
effective temperature of the system was constant, but 
the frequency would change hw e — hto — p. This effect 
is a result of our definition of the system-bath interac- 
tions. Those interactions consist of absorption or emis- 
sion of the bath quanta, and thus do not conserve the 
number of the bath particles. The Boltzmann exponent 
([37]) must account for energy gain or loss due to bath 
particles creation/annihilation. The remaining question 
is how the abstract operator algebra is represented in 
the Fock-Hilbert space of our multiparticle system. To 
answer this question we first introduce the operator 



D = B — A 2 /VN, 



(38) 



which commutes with A, A* , B and B^ . We then observe 
that (see appendix A): 
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• For each value of I — 0, or ! = 2,3,... there 
exist tin(1), so called, vacuum states \0i sv ) that 
are annihilated by the "jump" operator A — 
Er=o Vn + lglgn+i, 



A\0i sv ) = 0, 



(39) 



where the index I indicates the bare energy of the 
corresponding states (in units of fiu), 



22 nc4a„|0;, s ) = l\0i, s ), 

n=0 



(40) 



The state corresponding to I = contains all N 
atoms in the ground state. The states with higher 
energy are constructed as linear combinations of 
degenerated Fock states. There is only one state of 
I = 1 (N — 1 atoms in the ground state, and one 
atom in the first excited state), which is not anni- 
hilated by A; that is why there is no vacuum state 
with 1 = 1. The index s measures the number of ex- 
citations of in the state |0z s „) oc (D^) s \0i_2s0v) , 
and runs from to E(l/2), i.e. integer part of 1/2. 
Each of the states \0i sv ) is an eigenstate of D^D 
with an eigenvalue (4ls - 4s 2 - 2s + 2Ns)/N. The 
index v enumerates finally the states of energy I — 2s 
annihilated by D. Denoting by mjv(i) the number 
of states of energy I annihilated by D, i.e. the states 



D\0 



I0v) 



0. 



we obtain 



E(l/2) 

tin(1) — J]] iunQ — 2s). 

s=0 



(41) 



(42) 



The construction of the vacuum states is descibed 
in the Appendix A. Some other examples of the 
vacua are constructed in Appendix B. Each of the 
vacuum states is a linear combination of the acci- 
dentally degenerated energy eigenstates. The num- 
ber of accidentally degenerated states of the energy 
I in the N atom system, pn(1), is given by a solu- 
tion of the partition problem of the number the- 
ory Q] and is extravagantly large (c.f. Pn(1) — 
0(exp(7ri/2Z/3)) for I < N). The number of the 
vacua is given by njy(l) = Pn(I) — Pn(I — !)■ The 
very existence of these multiple vacuum states is 
thus a direct consequence and, at the same time, a 
signature of the accidental degeneracy |42| . 

• The vacuum states are orthonormal, (0i sv \0i' s i v ') = 
Su'6 ss '6 vv ' . 

• The Fock-Hilbert space of the system splits into an 
infinite number of Fock subspaces corresponding to 
each of the vacuum states. The Fock states in the 
(I, s, v)-th subspace are constructed as 



\kisv) 



\0lsv), 



(43) 



with k = 0, 1, They are also mutually orthonor- 
mal, and are also eigenstates of the energy operator 
with the corresponding eigenvalue (l + k). They are 
also highly degenerated (for k + I = k' + I'). In the 
following we shall use the notation w = [I, s, v). 

The dynamics exhibits in the Lamb-Dicke limit mul- 
tiple timescales. On the fastest time scale of the order 
of ry~ 2 it is nonergodic, i.e. it does not mix the different 
w-subspaces. After a short time of the order of the in- 
verse of T| all coherences between the \k w ) and \k' w ,) 
for k 7^ k' vanish. Within each w-subspace the system 
approaches the "thermal" equilibrium characterized by 
the density matrix diagonal in k, with zero off-diagonal 
elements, and the inverse temperature j3 e . The dynam- 
ics, however, cannot be reduced to a Pois son jump process, 
i.e. a sequence of random jumps between the various | k w ) 
states with the transition probabilities governed by the 
detailed balance conditions characteristic for the thermal 
equilibrium. The reason is that coherences corresponding 
to w 7^ w' , but k — kl do not vanish. 

The quantum mechanical density matrix in this (for- 
mally stationary) limit becomes a sum of diagonal canon- 
ical ensembles in each of the subspaces, accompanied by 
a sum of non-diagonal terms connecting different w and 
w 1 for the same fc's, 



= [1 - exp(-[3 e fiuj)}y^y^n w \k w )(k w \e- 



-f3 e hujk 



w k 



+ [l-exp(-/3 e W)] Yl Y. r ^'\ k ^^'\ e ~ ' hUJk - ( 44 ) 

w^w' k 

where the coefficients describe the populations of the cor- 
responding subspaces, and are defined as 



n w = Tr (P w p(0)P w ) . 



and the cumulative coherences, 



= ^<*gp(o)|Av)- 



(45) 



(46) 



Evidently, p exhibits nonergodic effects and depends ex- 
plicitely on the initial density operator. In the above 
expression P w = J2k l^tu)(^u»l denotes a projection op- 
erator onto the w-th subspace. Note that the definition 
(|H) implies that J2 W n w = 1. 

We stress that in the Fock basis spanned by the states 
\k w ) the matrix p is, in principle, not diagonal; moreover, 
in general, depending on the initial condition, it does 
contain coherences when represented in the Fock-Hilbert 
space corresponding to noninteracting atoms. That is the 
reason why the master equation in the latter basis can- 
not be reduced to the diagonal form in the Lamb-Dicke 
limit. That is why in order to arrive at such reduction 
additional assumptions have to be evoked that lift up the 
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accidental degeneracy, as anharmonicity of energy levels 
caused by anharmonicity of the trap potential or inter- 
atomic interactions. 



and the cumulative coherences 

r W w'(t) = y^(k w \p(t)\k v 



(49) 



VI. RESTORATION OF ERGODICITY 

From the previous Section it is clear that the break- 
down of ergodicity has only an approximate character 
since it is related to the lowest order dynamics in the 
LD expansion. It is natural to expect that Eq. (|44|) 
describes only a quasi-stationatry solution which is in- 
deed reached on a time scale l/r^.l/Tj ~ 0(oc 1/rj 2 ), 
but undergoes further slow evolution on a time scale 
l/r| 1} , l/r[ 1} ~ 0(oc l/r/ 4 ). On this slower time scale the 
ergodicity should be (at least partially) restored, and the 
transitions between the different w-th subspaces should 
become possible. As we shall see below, that is indeed 
the case. 

To this aim we consider the higher order corrections to 
the master equation (|2^) and treat them perturbatively 
within the standard adiabatic elimination method p3[ . 
We introduce the projection operator 

Vp{t) = n w (t) ]T \k w )(k w \e-P° n » k (l ~ e-^) 



w^w' k 



-j3 e hujk / 



\l-e- ° hw ), (47) 



where the populations of the corresponding subspaces 

n w (t)=Tr(P w p(t)P w ), (48) 



become now slowly varying functions of time. 

Introducing the complementary projector Q — 1 — V 
obtain 

Vp= V (4 X) + 4° ) V P + V ( C i ] + 4 1} ) Qp, (50) 
Q P = Q£ (0) Qp+Q + 4°) "Pp + ... (51) 

In the latter Eq. ( |5l| ) we have already employed the fact 
that rW/r ~ 0(cx rf), and neglected the higher order 
terms in rj 2 . Moreover, adiabatic elimination of Qp from 
Eq. (|5l]) introduces corrections of the order oc rf to Eq. 
( |50| ) for V p. Thus it may also be neglected in comparison 
to the leading term on the RHS of Eq. (|o|) . 

We obtain thus a relatively simple master equation 



(i) 



(52) 



From this master equation, after elementary algebra we 
obtain a set of closed rate equations for the popula- 
tions n w of the u>-th subspaces, and the coherences r ww i. 
These equations can be enormously simplified using the 
properties of the operators A, B, C, and their hermitian 
conjugates that are discussed in Appendix A. Amazingly, 
it is only the term C% which contributes in this order 
of the LD expansion; moreover, the equations for pop- 
ulations and coherences decouple (see Appendix C for 
details), 



n w =2T{ 1) (l-e-^)Y / E \(k w \B\k' w ,)\ 2 e-^ k ' n w , - (k w \B^ B\k w )e~^ k n 

k \k',w' 

+ 2T^(l - e~^)J2 ( E \(k w \B^\k' w ,\ 2 e-^ k 'n w , ~ (k w \BB^\k w ) e -^ k n 



k \k'w' 



2rW(l-e-^-)^ £ (k w \B\k' Wi )(k' W2 \B^\k w ,)e 



k \k',Wi,W2 



— \_(k w \B^ B\k Wl )v WlW / + (k Wl \B^ B\k w ')r WWl ^ 



-{3 e fiujk 



(53) 



2r W(l- e -^)^ £ (k w \B^\k' Wl ){k' W2 \B\k w ,)e 



-{3 e hu)k 



k \k',wi,W2 



^^[(KIBB^K,^^ + (k wl \BB^\k w ,) 



-/3 e huik 



(54) 







The above equations are further reduced inserting the 
unity between B and B\ and using the fact that A, and 
by definition do not couple the different w-subspaces, 
whereas B (B>) has the only relevant matrix elements 
between the different Zsw-subspaces for the (I ± 2, s ± 1, 
v-th.) and Isv-th subspace, for I = 2, 3, . . . (I = 0, 2, 3, . . .), 
s = 0, . . . , E(l/2), and k = k' = 0, 1, . . . (see Appendix 
C). 

From the above considerations we obtain the final form 
of the equations 



2ri 1} [fl 



+2 n w+ 2 



f w n w\ 



(55) 



2r 



(i) 



fw+2fw'+2 Tw-\-2w' ' -\-2 



2T 



(1) 



f iu fw ' — 2w ' - 



1 

2 ~2 



2 {fw + fw') r u 
(fw+2 + fw'+2) r u 



(56) 



where we have denoted w — (Isv), w±2 — (Z±2, s±l, v), 
with f w — (0 W _2\B\0 W ) . Note that f w can be assumed to 
be real and positive without loss of generality. As we see, 
both the populations n w , and the coherences r ww > fulfill 
a closed set of rate equations. The explicit expressions 
for the matrix elements involved in the above formulae 
are derived in Appendix C. 

The above equations provide the basic result of this 
paper. It shows that due to the presence of accidental 
degeneracy in the LDL the dynamics occurs essentially 
on several time scales. On the fastest time scale (gov- 
erned by C^) the dynamics is nonergodic and consist in 
approach toward the canonical equilibrium states in each 
of the Z-subspaces, accompanied by creation of quasi- 
equilibrium coherences between the states corresponding 
to the same fc's, but different w and w' . On this scale 
the populations of each of the w-subspace, as well as the 
cumulative coherences r ww i may be regarded as constant, 
and therefore the values of these coherences, as well as 
the populations in each of the subspaces depend on the 
initial conditions. In the higher order of expansion (on 
the time scale 1/r] 2 times longer), the mixing between dif- 
ferent w-subspaces becomes possible. This mixing leads 
to a partial restoration of ergodicity. In the example con- 
sidered here this restoration is not full, however, because 
as easily seen from Eq. ( p5| ) the w-subspaces with differ- 
ent v still do not mix. The reader can easily check that 
the further mixing of the odd and even subspaces does 
take place in the next order of the Lamb-Dicke expansion 
(for instance due to term containing bi-products of the 
operators C, and C 1 *). 

The stationary state that is reached on the slower time 
scale is easy to find since Eqs. ( pS] ) fulfill the detailed bal- 
ance conditions, whereas the cumulative coherences are 
damped to zero, as demonstrated in Appendix D. We 
obtain that 



so that 



with 



n2l+m,l,v — e «mO»i 



e = /3- Pn/2Kw. 



(58) 



(59) 



The ratio of n m o Vl and n m >ov' remain undetermined in 
this order. Note that the reason why [3' e ^ [3 e is that 
both temperatures correspond to the processes that in- 
volve single bath quantum absorbtion or emission, but 
different energy changes (by ftio, or 2huj, respectively). 
Indeed, it is elementary to check that the stationary di- 
agonal matrix elements of the desity matrix are propor- 
tional to the corresponding Boltzmann factors, 

{k2i+m,l,v\p\k2l+m,l,v) oc e'^ k ~ 2 ^ l n m0v 

= <fc 2 / +m , ; ,„|e-' 3 =^^- m+/3 ^ tA/2 n m o.|A ; 2 ;+m ,/^>- (60) 

Finally, it is easy to check by substitution in the Liou- 
villian (|l4|) that for /i = the steady state to all orders 
in the LD expansion is precisely p cx e~@ HA , which is di- 
agonal in the original basis. Obviously, the steady state 
solutions obtained in the first, and the second order of 
our expansion are compatible with such a steady state. 



VII. CONCLUSIONS 



ni+2,s+i,v 



ze riisy. 



(57) 



In a series of papers ]2 2| |24| . |2 7| | we have studied in detail 
the quantum dynamics of bosonic and fermionic gases of 
cold atoms in traps in the absence of accidental degener- 
acy. We studied various cooling mechanism, and various 
limiting cases. In this paper we have presented the so- 
lution of the corresponding problem accounting for acci- 
dental degeneracy effects. We have studied interactions 
of a gas of trapped atoms with a heat bath in the Lamb- 
Dicke limit using the master equation approach. We have 
demonstrated that the system approaches an equilibrium 
on two (or more) distinguished time scales, and that the 
dynamics has the corresponding number of stages. At 
each stage a quasi- equilibrium state within appropriately 
determined subspaces of the Hilbert space is reached. 
This quasi-equilibrium corresponds to a canonical ensem- 
ble resticted to the appropriate subspace, and character- 
ized by an effective temperature determined by the tem- 
perature and the chemical potential of the heat bath. In 
the next stage thcrmalization between the groups of sub- 
spaces occurs leading to another quasi-equilibrium in the 
larger subspaces, and so on. 

We would like to stress once more that the problem 
considered in this paper is quite general. Atomic traps, 
although frequently anisotropic (see for example Ref. Q), 
can be designed to be harmonic with a very good accu- 
racy. A cooled atomic sample in such a harmonic mi- 
crotrap will necessarily exhibit the effects of accidental 
degeneracy regardless of the method used for its coolingl 



10 



One may question the generality of our results, since 
we have used the Lamb-Dicke expansion, and at the same 
time neglected in this paper atom-atom collisions, as well 
as atom-atom interactions mediated by the coupling with 
the bath. Such processes (described by the Hamiltonian 
H' A _ A) see Section III) lead evidently to shifts of the 
atomic energy levels, and will, in principle, lift up the 
accidental degeneracies. As we argued in Refs. [f22|~p4|, 
as long as the number of atoms in the microtrap is not 
too large, those shifts remain small and can be treated 
perturbatively. The system will then still exhibit the ef- 
fects of accidental quasi-degeneracy. We stress that the 
theory developed in this article is valid for arbitrary num- 
bers of atoms, and in particular it is for two atoms. Using 
far-off-resonance dipole traps [Q , or loading atoms to a 
single minimum of a dark optical lattice (see Ref. |(|(b)) 
it should be accesible experimentaly to confine several 
atoms in the trap of the size a ~ 0.1 — 0.05/xm. That 
implies validity of the LD expansion for the laser colling 
case (see (p2j). Similarly, one can use a cooled atomic 
gas close to, or below the condensation point as a bath 
in the symphatetic cooling case. In the conditions of the 
experiments of Refs. S and [El that implies de Broglie 
wavelength of the order of /im, and thus validates the LD 
expansion. Using Bogoliubov-Hartree theory it is possi- 
ble to estimate perturbatively that the energy level split- 
tings in a " band" of the quasi-degenerated states due to 
atom-atom collisions will be in such a case of the order 
of N(a sc /a)/ \/Nd, where a sc is the scattering length of 
the system atoms, and Njj is the number of levels in the 
band. Note that Nd increases dramatically for higher 
excited levels. We see with a sc = 5 nm, our theory 
should work for N up to ~ 20 even in the worst case 
when Nd = 2. Note that the cooling of 20 atoms to 
the ground state of a harmonic trap might be a very in- 
teresting task for the rapidly-developing research field of 
quantum information processing. Additionaly, we want 
to recall at this point that as pointed above, the external 
modification of the s-wave scattering length via Feshbach 
resonances has been demonstrated, been experimentally 
feasible the achievement of a quasi-ideal gas. 

The main physical results of the paper are thus the 
following. We have been able to treat analytically the 
quantum dynamics of an ideal gas of N atom in the LD 
limit. We have shown that the dynamics naturally splits 
into two parts: a fast part, during which coherences be- 
tween the degenerated states are preserved, and a slower 
part, during which thermal equilibrium is achieved. Even 
though, the ideal case considered is not realistic (at least 
without external modification of the scattering length), 
we think i) that it provides a lot of insight into more 
realistic situations; ii) it is, to our knowledge, one of 
the extremely rare examples of soluble quantum dynam- 
ical problems in the area of statistical physics. The 
method that we developed, and results can be carried 
over to more realistic situations concerning cooling of 
small atomic samples (~20 atoms) in microtraps. Such 
situations are not far from the reach of present exper- 



iments. The calculations for such a case should follow 
exactly the lines described in this paper, with the only 
difference that the parts of the Liouville equations de- 
scribing the atom-atom interactions that lift up the ex- 
act degeneracies should be included into the corrections 
to the £(°) . In another words, they should be treated just 
like the corrections to the lowest order term in the LD 
expansion have been treated in this paper. It is obvious, 
that as long as the splittings of the quasi-degenerated 
levels will remain small relative to Tilo, such realistic sys- 
tem will exhibit basic effects presented in this paper, i.e. 
step- wise dynamics on the two time scales. 

The quantum dynamics of samples of cold atoms ex- 
hibits, in our opinion, an enormous reachness of inter- 
esting physical and mathematical phenomena, such as 
multistable, exotic stationary states, multistage dynam- 
ics etc. The present paper is another example of this 
reachness. Further studies are, however, required to get 
more understanding of this new physics, including for 
instance developement of other statistical physics tools 
p8| (c), such as diffusion equations, hydrodynamic limits 
etc. 

After this work was finished we have learned from T. 
Fischer, K. Vogel and W. Schleich that similar algebra to 
the one used by us appears in the problem of the cool- 
ing of a sample of bosons with a simple particle reser- 
voir Ga|. We thank Yvan Castin, Jean Dalibard, Ralph 
Dum, T. Fischer, K. Vogel, and P. Zollcr for enlighting 
discussions. We acknowledge the support of the Deutsche 
Forschungsgemeinschaft (SFB 407), ESF PESC Proposal 
BEC2000+, and TMR ERBXTCT96-0002. 



APPENDIX A: THE STRUCTURE OF THE 
FOCK-HILBERT SPACE 

The matrix elements of the operators B, and C can be 
calculated using elegant algebraic methods. To this aim 
we first observe that the operators in question fulfill the 
commutation relations 



[A,B] = 0, 




(Ala) 


[A,C} = B/(2v 


/ iV), 


(Alb) 


[A,Bl 


= 2A+/V 


% 


(Ale) 


[A,cr 


= E/N H 


-1/2, 


(Aid) 


~B,B r 


= AE/N 


+ 2, 


(Ale) 


B,E 


= 2B, 




(Alf) 


C,E 


= c, 




(Alg) 


A, E 


= A, 




(Alh) 



with E = X)^Lo na n a n denoting the normalized energy 
operator. 

It proves to be very useful to introduce the operator 
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D = B — A 2 /^fN. 
This operator fulfills 



(A2) 



[D,A] = 0, (A3a) 
[D, A f ] = 0, (A3b) 
[D,D^\ = A(E- A A)/N + 2(1- 1/N), (A3c) 

Since the operators A' A, and D'D commute, it is useful 
to characterize the multiple vacua in terms of the eigen- 
values of these two hermitian operators. 
Let \0iov) denote the states that fulfill 



A|0 i0 „) = 0, 
D\0 lOv ) = 0. 



(A4a) 
(A4b) 



There are mjv(i) such states, and the index v enumerates 
them. Note that the states 

\0 lsv ) = (i>t)- j 0l _ a , j0 ^>/| | (2?+)-|0,_3», 0l «> 1 1 , (A5) 

have energy I, are annihilated by A, and are eigenstates 
of D^D with the eigenvalue 



s-l 



W - 2s + 2s') /N + 2(1 - 1/N)} . 



(A6) 



s'=0 



In the subsequent energy sectors we have thus the vacua: 
|0ooi), |0 2U >, |0soi), |0 42 i), IO401), IO511), |0 50 i), IObsi), 
|0 6 n), |0 60 i), |0 602 ), etc. 
The Fock-Hilbert space is spanned by the vectors 



(A^) k 

\h sv ) = [ -^r\o ls 



(A7) 



APPENDIX B: CONSTRUCTION OF VACUUM 
STATES 

We have seen in Appendix A that the vacuum states 
can be constructed by applying the operator D* consec- 
utively to the states |Oz 5 o,w) - I n this Appendix we present 
explicit construction of another family of the vacuum 
states that are annihilated by the operator A. In fact 
we consider a more general case with 



A = ^2 A n ala n+1 . 



(Bl) 



71=0 



Such defined operator reduces to the one given by Eq. 



(27a) if we put A n = \Jn + 1. The vaccum states fulfill 
A[0i) = (B2) 



There is one obvious solution of the above equation 
which describes the global ground state 



Apart from that, for I = 2, 3, ... we define 



l-l 



|0l)= Y>a l m al +1 _ m \N-m,m- 1,0, ■■■) 

TTL—1 

+a\\N -1,1,0,...). 



(B4) 



From the above definition it is clear, that different vacua 
are orthogonal. Applying A to the above expression after 
elementary algebra we derive the recurrence formulas for 
the coefficients 



A 



l — in 



valid for m = 1, 1 — 2, and 



l m+l 



(B5) 



H-l 



A VWN -l + l , 



(B6) 



From the above expression it is easy to construct ex- 
plicitely the corresponding vacuum states. The value of 
a\ can be conveniently chosen for normalisation of the 
states. 



In the case considered in this paper (A n — \Jn + 1 ) 
the first few normalized vacuum states are: 



|0 2 ) = -=(|A-2,2,0,. 



VA^T| N -1,0,1,0,...)) , 
v/8 



(B7) 



VA 2 + 3A-2 
V3(JV-2) 



(|A-3,3,0,...) 



N -2,1,1,0,. 



V(N-l)(N-2) 



2V2 



I N- 1,0,0,1,0,...) 



(B8) 



|0 4 ) 



y/N 3 - 5N 2 -3N+ 21 

j2(N-3), 

iV-3,2,1,0,. 



(|A-4,4,0,...) 



2y/(N-2)(N-3) 



|A-2,1,0,1,0,. 



y/(N-l)(N-2)(N-3) 



\N- 1,0,0,0,1,0,...) (B9) 



|0o) = |JV,0,0, 



(B3) 



etc. 
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APPENDIX C: CALCULATION OF THE 
MATRIX ELEMENTS 

Let us first consider the operator B, and derive the 
explicit expressions for the martix elements 



flsv — \0l-2s-lv\B\0l. 



(0i- 2 .-it,|-D|0i. 



that enter Eq. (|55|). Note that the coefficients fi sv are 
real, since the matrix elements of the operator A in the 
non-interacting atom basis are real (see Appendix A). 
Moreover, without any loss of generality /; s „'s may be 
assumed to be non-negative. In the following I will skip 
the index v which is not affected by the dynamics. 

Since we know that B\0i s ) cx |0;_2s-i) for I > 2s, we 
can write 



B\0 ls ) = / Ja |0,_ a .-i) 



(C2) 



Similarly, using the commutation relation ( |Alq ) we can 
write 



B f |0 



!— 2s— lv, 



= fis\0is) + V2/N\ 



fl = (0,_ a , B _i|SBt|0,_ 2 ,,_i) - 2/N. 



(C3) 



(C4) 



From the above expressions using the commutation re- 
lation (Ale) we obtain 



./z+2, + i = fl + 41/JV + 2 - 2/N. (C5) 

The above recurrence can be easily solved yielding 

ff= (2 - 2/N)s + 4((J - 2s)s + s(s - 1))N. (C6) 

since /;-2s,o = 0. 

In general, we may write 



B\k ls ) = f u \h-2,.-i) + VHk-l]/N\(k - 2) u ). (C7) 

The above formulae provide a very efficient method of 
calculating all of the matrix elements of the operators B 
and . 

It is a little more tedious to derive corresponding for- 
mulae for the operator C. From Eq. ( Ala| ) we obtain 



AC\0 ls ) = 



1 



2VN 



/JO 



Z-2,s-l/ 



(C8) 



so that 



E(Z-l/2) 

C|0 is ) = - 7 =/ is |l / _ 2 , s _ 1 )+ w lA®i-i,s>)- (C9) 



2ViV 

The coefficients wi s > — (0;_i^_i|C|0/ s ) may be also re- 
garded to be real, and can be, for instance, determined 
directly from the definitions of the vacuum states. 
In general, we can write 



C\k ls ) = ^^fulk + h^i 



2VN 

E (7-1/2) 

+ n^\ki-i,^)+vi,(k)\k-ii.), (Cio) 



(Cl) with 



v ls (k + l) = 



-vis{k) 



1 



Vfc+T 



N 



I 1 
- + 2 



(Cll) 



and vi s (0) = 0. The above formulae allow for very effi- 
cient calculations of the matrix element of the operators 
C and C', provided the states |O m ou) are known. 



The expression (CIO) implies immediately that £^ 
does not contribute at all to the final equations ( |54| ) 
and (^4|). Since the matrix p is diagonal in the k in- 
dex, whereas the operators A and A^ change k to k — 1, 
and k + 1 respectively, only those parts of the operators 
C and that change k ± 1 back to A; could contribute. 
It is evident from Eq. ( plO ), however, that these parts 
of C and do not change Isv. Therefore, their contri- 
butions to Eqs. ([34]) and ( j54| ) vanish identically, due to 
the trace-like sums over A: appearing on the right hand 
side. 

Similar considerations show that there is no mixing of 
the populations and the cumulative coherences in Eqs. 
( f54| ) and (p4]). Let us, for example, consider Eq. (H) for 
the populations ni sv . As in the previous case, the con- 
tributions of the parts of the operators B and B' that 
do not change Isv vanish identically, due to the trace- 
like sums over k appearing on the right hand side. The 
parts of B and B^ that change I by =f2, and s by =f1 
do contribute, but they can only transform the parts 
of the density matrix proportional to \ki sv ){ki> s i v >\ into 



x 



\h sv ){k v 



or oc K;±2, s ±l 



,v){kl 



±2,s'±1,d' 



they 



can only lead to couplings between the populations ni sv 
and ni±2,s±\,vi or between the cumulative coherences 
Tisv,i' s 'v' and ri±2 :S ±i,v,i'±2,s'±i,v' ■ 



APPENDIX D: DECAY OF COHERENCES 



In this Appendix we keep a single index w — (l,s,v), 
and denote w ± 2 = (I ± 2, s ± 1, v). In order to prove 
that the cumulative coherences decay to zero on the slow 
time scale, we rewrite Eq. (p6|) in the form 



■DB 



■NEG 



(Dl) 



where the first term 



T'ww' — ' [f w +2.fw'+2l r w+2w'+2 — fwfw'fww'] 

"H2T 1 [fwfw'Tw — 2w '-2 ^ fw+2fw'+2 'WW' } , (D2) 
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corresponds to a set of kinetic equations with (positive) 
rates that fulfill detailed balance conditions. The matrix 
that enters the right hand side and generates the evolu- 
tion is therefore evidently non-positively defined, and has 
exactly one eigenvector with zero eigenvalue, correspond- 
ing to the stationary solution of the Boltzmann-Gibbs 
form. 



The second term in Eq. (Dl) 



■NEG _ 

' 1/11!]' I 



l 



{fw fw' fwfw') r u 



-2T 



(1) 



2 (fw+2 + fw'+2 — fw+2fw'+2) r u 



(D3) 



corresponds to a set of simple decay equations, with the 
rates which are evidently positive, since |(/^ + f^i) — 
fwfw' is strictly greater than zero for w 7^ w' . 

The full dynamics of the cumulative coherences is thus 
generated by the sum of the two martices, one of which is 
non-positively defined, and the other being strictly neg- 
atively defined. The sum itself must therefore be nega- 
tively defined, ergo it generates the decay to zero. 
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